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■ Cold atom systems held in optical lattice potentials are 
proving to be a remarkably versatile means of investigating 
^ coherent quantum phenomena. A particularly interesting 
Ol"ecent development has been the use of these systems to 
study synthetic gauge fields. In Ref. [1] a scheme was 

S proposed to simulate a uniform magnetic field by applying 
a periodic driving to the optical lattice, and employing the 
effect known as "photon assisted tunneling" to generate 
£j the required Peierls phases. Here, however, we show that 
Othe analysis of this system was incorrect; contrary to the 
i ^ c entral claim of the paper, a driving of the kind described 

cannot produce a uniform synthetic magnetic field. 
C^l , For simplicity, we begin by considering a particle hop- 
^ bing on a one-dimensional lattice, described by a tight- 
qq binding model 

ffo = -^(|m)<m+l| + tf.c.) > 

cn : 

where \m) arc Wannier states localized on sites to, and 
<7 is the hopping between nearest neighbors. Oscillating 
. .the position of the optical lattice, x — > x + xq cos(wt + <f>), 
^"produces an inertial force in the rest frame of the lattice 
Fj = — Mxoui 2 cos(wi + </>), which can be described in 
jjj terms of a potential function 
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V(x,t) 



(2) 



This allows the driving to be included in the Hamiltonian 

as 



H = H Q + V(x,t) 

= H + aF^ cos(uit + cj>) \m)r 



(3) 



where a is the lattice spacing and F u = Mxqlo 2 . In addi- 
tion to the oscillation the lattice can also be subjected to 



a constant acceleration [5] , which adds a constant compo- 
nent F to the total lattice potential 

rn 

+ a(F + F u cos(cjt + 4>)) \m)m(m\ . (4) 

m 

It can be readily seen that this Hamiltonian mimics that 
of a charged particle subjected to a uniform electric field, 
with F and F^ playing the roles of the DC and AC com- 
ponents respectively. 

When the static component matches the driving fre- 
quency, F — us, a straightforward Floquct analysis 3,4 
reveals that in the high-frequency limit the dynamics of 
the driven system can be described by an effective static 
Hamiltonian 
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— (\ m + l)He # + H.C. 



JJi(F u /uj) 



(5) 



where J\{x) is the Bessel function of the first kind. The 
presence of the phase factors, e ±l< ^, in the effective Hamil- 
tonian provides the motivation for attempting to use this 
kind of driving to generate synthetic gauge potentials. 

A two-dimensional square lattice threaded by a uniform 
magnetic field can be described in the Landau gauge A = 
4(0, x,0) by 
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^2{\l + l,m)(l,m\ + H.c.) 



(6) 
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Fig. 1: Evolution of the initial state ip(x,y) = 
N (exp[-2(a; 2 + y 2 )] + exp[-2(a; 2 - (y - 8) 2 )]) . (a) Den- 
sity distribution of the system at t = 30T after evolving in 
the presence of a uniform magnetic field (Eq(SJ. (b) As in (a), 
but for a periodically-driven system. The peak centered on 
the origin shows a similar evolution to the exact case, but the 
other peak is instead compressed into a narrow strip parallel 
to the j/-axis, due to the renormalization of J x . Parameters 
are u = 19.9, F/u = 1, F u /ui = 0.2, and 4> = tt/3. 



where J x /J y are the tunneling amplitudes in the x/y di- 
rections, and (I, m) are the coordinates of the lattice sites. 
In Ref. [I] it was proposed that this pattern of Peierls 
phases could be mimicked by adding the periodic driving 
potential 



V{x, y,t) = y(F + F u cos(wt - x</>)) 



(7) 



to the two-dimensional tight-binding model. With this 
driving potential, the force experienced in the rest frame 
of the lattice (F = — VV^) is given in the continuum limit 
by 



F, 
F„ 



-d x V 
-d y V 



[c/>F u sin (ut - x<f>)] y (8) 
[F + F u cos (ojt - x4>)} . (9) 



Comparing ([§]) with Eqs0] and [5] indicates that the tun- 
neling in the y-direction will indeed acquire the correct 
phases to reproduce the Landau gauge. However, it is 
also apparent that tunneling in the x-direction will ex- 
perience a periodic driving (|5]) too, and so will also be 
renormalizcd. This may appear unexpected, as this renor- 
malization is occurring for tunneling perpendicular to the 
shaking of the lattice, but is a consequence of the explicit 
x-dependence of V(x, y, t). This effect was not considered 
in Ref. [T|, as was implicitly assumed to be equal to 
zero. Consequently the results obtained in Ref. [I] are 
only valid for F x ~ 0, that is, in a sufficiently small region 
around y = 0. 

Performing a Floquet analysis reveals that while J y is 
renormalizcd as in EqJSJ the the effective rr-hopping is 
given by 



jf = J x Jo sin [0/2] — 



(10) 



this effect. The system is initialized in a superposition of 
two narrow Gaussians, one centered on the origin, and the 
other on (0,8). FigJT^, shows the system evolving under 
the action of Eq[6]for a magnetic flux of (f> = 7r/3, while 
in FiglT}} we show the result of evolving the system under 
the time-dependent driving potential ([7]) . We can see that 
the peak centered on the origin behaves similarly in both 
cases, indicating that near this point the driving indeed 
mimics a magnetic field with reasonable accuracy. The 
behavior of the other peak, however, is strikingly differ- 
ent. The suppression of tunneling produced by the Bessel 
function squeezes the wavepacket into a narrow band, with 
a very different dynamics to that of the true magnetic flux. 

In summary we have shown that the scheme proposed in 
Ref. pP cannot produce a uniform magnetic field. At best 
it can produce a field that is approximately constant over 
a finite number of lattice spacings, which severely limits 
its applications. We have demonstrated this by an explicit 
calculation of the effect of the driving potential. This re- 
sult is also consistent with the basic observation that a 
vector field with a non-zero curl cannot be described by 
a conservative potential. More complicated experimen- 
tal setups are required to avoid this drawback; for exam- 
ple introducing additional optical lattice potentials [5] to 
compensate for the transverse renormalization (|10j) . 
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Thus away from y = 0, J x will be suppressed by the ze- 
roth Bessel function. In FigQ]we show the magnitude of 
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